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Abstract 

The extracted value of the relic density has reached the few per-cent level precision. One can there- 
fore no longer content oneself with calculations of this observable where the annihilation processes are 
computed at tree-level, especially in supersymmetry where radiative corrections are usually large. Im- 
plementing full one-loop corrections to all annihilation processes that would be needed in a scan over 
parameters is a daunting task. On the other hand one may ask whether the bulk of the corrections are 
taken into account through effective couplings of the neutralino that improve the tree-level calculation and 
would be easy to implement. We address this issue by concentrating in this first study on the neutralino 
coupling to i) fermions and sfermions and ii) Z . After constructing the effective couplings we compare 
their efficiency compared to the full one-loop calculation and comment on the failures and success of the 
approach. As a bonus we point out that large non decoupling effects of heavy sfermions could in principle 
be measured in the annihilation process, a point of interest in view of the latest limit on the squark 
masses from the LHC. We also comment on the scheme dependencies of the one-loop corrected results. 
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1 Introduction 



With barely lfb - of data, the LHC is pushing many hitherto popular, though naive, extensions 
of supersymmetry to the corners of high masses [lj while leaving some hope for a discovery of a 
rather light Higgs that could still be compatible with supersymmetry [2j. Before this very recent 
paradigm, susper symmetric models (and most models of new physics for that matter) were very 
strenuously constrained to a thin sliver in parameter space, most notably from the very precise 
measurement of the dark matter relic density that has now reached the few percent level and that 
will get even more precise in the future, hence cornering even further model building. Combining 
the results of the 7-year WMAP data [3] on the 6-parameter ACDM model, the baryon acoustic 
oscillations from SDSS[4J and the most recent determination of the Hubble constant [5] one[6j 
arrives at ^cdm^ 2 = 0.1123±0.0035, where ^cdm is the density of cold dark matter normalised 
to the critical density, and h is the Hubble constant in units of 100 km s Mpc -1 . One has 
reached a precision of 3%. The data from the LHC does not infer that dark matter within 
supersymmetry, exemplified most nicely through the neutralino, the lightest supersymmetric 
particle (LSP), is of order of ITeV or so, this is just a limit on the coloured constituents of the 
model. As for the Higgs, were it not for the large radiative corrections on the mass of the lightest 
member, supersymmetry would long be a forlorn construct. The Higgs is the most prominent 
example where radiative corrections are far from negligible in supersymmetry, yet practically all 
analyses that aim at constraining the parameter space of the MSSM through the relic density 
are based on tree-level cross sections of the annihilation processes entering the predictions of this 
quantity which as stressed is experimentally given within the percent accuracy. Only seldom do 
some analyses assign a theory uncertainty to these annihilation cross sections, an uncertainty due 
essentially to the fact that higher order loop effects are not known. This uncertainty, in the rare 
case where it is taken into account, is however assumed to be invariably the same whatever the 
nature of the dominant process and the composition of the LSP. The reason the loop corrections 
are ignored, irrespective of the model specified, is that the calculation of the relic density requires 
most often the evaluation of a large number of processes. Most analyses are done with public 
codes [3 [HJ |9] based on tree-level calculations. Computations of the relic density at one- loop have 
now been achieved for quite a few channels |10[ [Til \12\ [TBI Q3] and tools exist now to perform in 
principle any calculation of the relic density beyond tree-level amplitudes thanks to the recent 
development of adapted automation tools [15]. Beside the findings [TO] [TTJ [12j [131 E] that these 
corrections are important, the improvements have not percolated to most analyses. It must be 
said that these calculations do involve some non trivial issues about the renormalisation of the 
MSSM and more generally techniques for one-loop integrals that certainly require expertise. The 
other reason is that even when they could be implemented, they are still certainly extremely 
CPU time consuming, forbidding hence any attempts of fits, likelihood search, and in a more 
general context any sampling of the parameter space, especially if one takes into account the 
fact that the MSSM is more than liberal with unconstrained parameters. Yet, apart from 
providing a more precise prediction, one-loop calculations can probe higher masses, a situation 
akin to the precision electroweak observables and their sensitivity to the top and Higgs mass. 
For example, non decoupling effects termed in analogy with electroweak SM observables, super- 
oblique corrections have been revealed in one- loop calculations of supersymmetry observables [16[ 
1X71 118]. An example in view of the recent findings of the LHC is that super heavy squarks leave a 
non negligible imprint on many observables in particular the annihilation cross sections involving 
dark matter. 

The aim of this paper is two-fold. First, to stress again the importance of the loop corrections 
for the relic density and show again that even when a loop calculation is available, there still 
remains in some cases an uncertainty that pertains to the choice of the renormalisation scheme. 
The second and more detailed aspect is to discuss whether an approximation to the one-loop 
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calculation can be found. We aim at implementing a universal correction through effective 
couplings of the LSP and check its validity against a complete one-loop calculation. If such an 
approximation is possible and general enough it could be implemented in existing codes (based 
on tree-level cross sections) calculating the relic density. Such was the case with the inclusion of 
the Sommerfeld effect [19J in the case of coannihilation or processes dominated by Higgs exchange 
for which Arrib[20\ corrections are included [7 1 . 

This preliminary study takes a simple process, namely X1X1 ~~ * ^ + 1 1 ^ as a testing ground. 
The aim of this study is not to find a good scenario that returns the correct actual value of 
the relic density but to try to unravel some general common features of the loop calculation 
to improve the predictions of the relic density. The aim is to rather find out whether one can 
improve on the tree-level calculation by introducing effective couplings of the LSP that could 
be used for any process. As we will see, though at first sight naive, the process X1X1 ~^ 
embodies the three types of couplings of the neutralino: to ///, gauge bosons and Higgses. 
Here we cover the bino and the higgsino case. One might argue that the bino case corresponds 
to what was referred to as the bulk region in the constrained MSSM and is largely ruled out, 
whereas for a higgsino this would not be the dominant process. As we have just stressed the 
aim is not to strive to find a good scenario. Besides, it suffices to change the cosmological 
ingredients [21] entering the calculation of the relic density to revive the so called bulk region. 
We will see that by considering these few simple examples the conclusions about the efficacy 
of the effective coupling is quite different. Moreover, there is no need here to convert a full 
corrected cross section into a relic density value, we rather take for all the models we study the 
annihilation cross section for an energy corresponding to a relative velocity of v = 0.2, typical 
for the relic density calculation. As known, for zero relative velocity the process enjoys chirality 
suppression which is lifted at higher order through gauge boson radiation, however the effect on 
the relic is totally marginal|10]. 

The paper is organised as follows. We first briefly describe the ingredients necessary to 
perform a one-loop calculation in supersymmetry covering both automation, renormalisation 
and renormalisation schemes, it is in this section that we will write down the effective universal 
neutralino couplings as well as some definitions. Section 3 contains our main results. After a 
few definitions we first study the case of a bino-like neutralino before addressing the case of a 
higgsino-like LSP. We also quantify the possible uncertainties due to the scheme dependence. 
We conclude in Section 4 by some general observations. 

Throughout the paper we use some shorthand notation for angles. Generically C0,sg,tg stand 
for cos 9, sin 0, tan 8. The weak mixing angle dyy is defined as cos 6w = c\y = Myy/Mz where 
M]y is the W mass and Mz the Z mass. 

2 Calculations, renormalisation, schemes at full one-loop 
2.1 Tree- level considerations 

At tree-level, see Fig. [TJ XiXi ~^ ^ + ^~ proceeds through i) t-channel smuon exchange dominated 
by a jiR in the case of a bino-like since it has the largest hypercharge, ii) a Z exchange which, 
on the other hand, is suppressed for the bino Hi) Higgs exchange but this is small in view of 
the Yukawa coupling of the muon. Therefore, as advertised, all types of couplings for the LSP 
are present: to fermions in the x?W coupling, gauge bosons in the X1X1Z, and Higgs scalars 
such as XiXi-A (A is the pseudoscalar Higgs). It is through the choice of a hierarchy in the 
set Mi,M2,fi that we can largely define the nature of the LSP. Numericaly speaking we call a 
neutralino pure or almost pure when its mixing to the specified species is over 99%. 
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Figure 1: Tree-level diagrams contributing to X1X1 ~^ / u+ / u ■ (a) is the t-channel fl exchange, 
(b) is the Z exchange and (c) an example of a Higgs exchange. 

2.2 Renormalisation and loop corrections, general considerations and issues 

2.2.1 Set up of the automatic calculation: SloopS 

One-loop processes calculated via the diagrammatic Feynman approach involve a huge number of 
diagrams even for 2—7-2 reactions, especially in a theory like supersymmetry. Performing a full 
one-loop calculation by hand without automation is practically untractable. Our exact full one- 
loop calculation is done with the help of the automated code SloopS. SloopS is an automated 
code for one-loop calculations in supersymmetry. It is a combination of LanHEP [22j, the bundle 
FeynArts[23j, FormCalc[23] and an adapted version of LoopTools[251 [SB]. LanHEP deals with 
one of the main difficulties that has to be tackled for the automation of the implementation of 
the model file, which is entering the thousands of vertices that define the Feynman rules. On the 
theory side a proper renormalisation scheme needs to be set up, which then means extending 
many of these rules to include counter-terms. This part is done through LanHEP which allows 
to shift fields and parameters and thus generates counterterms most efficiently. The ghost 
Lagrangian is derived directly from the BRST transformations. The loop libraries used in 
SloopS are based on LoopTools with the addition of quite a few routines in particular those for 
dealing with small Gram determinants that appear in our case at small relative velocities of the 
annihilating dark matter, and even more so of relevance for indirect detection |26j. 

2.2.2 Renormalisation 

In SloopS all sectors of the MSSM are implemented through a one-loop renormalisation. This is 
explained in details in |10|, [27 ], I28| . I13J . Here we only briefly sketch the renormalisation procedure. 
We have worked, as far as possible, within an on-shell scheme generalising what is done for the 
electroweak standard mo del [29]. 

i) The Standard Model parameters : the fermion masses as well as the mass of the W and Z 
are taken as input physical parameters. The electric charge is defined in the Thomson limit, 
see for example [29]. The light quarks (effective) masses are chosen [30J such as to reproduce 
the SM value of a -1 (M|) = 127.77. This should be kept in mind since one would be tempted 
to use a DR scheme for a, defined as Mz, to take into account the fact that dark matter is 
annihilating at roughly the electroweak scale, so that a((2m^o) 2 ) is a more appropriate choice. 
One should remember that the use of a((2m^o) 2 ) instead of the on-shell value in the Thomson 
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limit would correct the tree-level cross section for XiXi ~~ ► by about 14%. As we will see 

and have reported somewhere else for other processes this running does not, most of the time, 
take into account the bulk of the radiative corrections that we report here. Therefore for further 
reference, let us introduce the correction due to the running of a, 

A a = ^L^£ = 2Aa, (1) 

where the cross section oq is the tree-level calculated with ao = a(0) = 1/137.0359895 whereas 
o" ae ff is the tree-level with ao -4 a e ff. (Q 2 ) = ot{Q 2 ) = olq/{1 + A.a(Q 2 )),Q = 2m^o. With our 
input parameters Aa(M|) = 0.06. In the running we allow for all charged particles including 
the W boson contribution, the top and the sfermions and the charginos, though for the light 
LSP scenario we consider these added contributions are very smalQ 

ii) The Higgs sector : The pseudoscalar Higgs mass M\ is used as an input parameter while 
insisting on vanishing tadpoles, tp, which at tree-level is the ratio of the two expectation values 
of the Higgs doublets, can be defined through several schemes: 

- In the DCPR scheme [3T|l32j tg is defined by requiring that the (renormalised) A°Z transition 

vanishes at Q 2 = M^ . 

- a DR definition where the ta counter-term, 5tp, is defined as a pure divergence leaving out 

all finite parts. 

- a process-dependent definition of this counter-term by extracting it from the decay ^4° — > 

t + t~ that we will refer to as A TT for short. This definition is a good choice for the gauge 
independence of the processes. 

- an on-shell definition with the help of the mass of the heavy CP Higgs H taken as input 

parameter called the MH scheme from now on. We have reported elsewhere that this 
scheme usually introduces large radiative corrections. 

These schemes are critically reviewed in [27] . By default we use the DCPR scheme but when 
quantifying the effect of the scheme dependence on ta we also use the DR and MH scheme. 

iii) The sfermion sector : For the process at hand only the smuons parameters require renormal- 
isation. For the slepton sector we use as input parameters masses of the two charged sleptons 
which in the case of no-mixing define the R-slepton soft breaking mass, M^ R and the SU(2) 
mass, Mfi L , giving a correction to the sneutrino mass at one-loop. Though not needed here, in 
the squark sector each generation needs three physical masses to constrain the breaking param- 
eter Mq l for the SU(2) part, M^ fl , Mj for the R-part. See[28] for details. 



iv) The chargino/neutralino sector. First of all, for the neutralinos at tree- level the physical fields 
X®, i = 1, • • • ,4 are obtained from the current eigenstates (tp n ) = (B°, W°, Hf, H®) through a 
unitary complex matrix N 

X° = (2) 

iV diagonalises the mass mixing matrix Y in the neutralino sector, see |28j for details and 
conventions. Although only Xi enters our calculations we do need to fix all the elements that 
define its composition and hence couplings. For this sector we implement an on-shell scheme 

x For the W boson contribution the self energy of the photon is calculated in a non-linear gauge[29] correspond- 
ing to the background gauge in order to maintain (7(1)qed gauge invariance. 
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by taking as input three masses in order to reconstruct the underlying parameters Mi, M 2 , [i. 
In SloopS [28] the default scheme is to choose two charginos masses m-± and m~± as input 
to define M 2 and /x and one neutralino mass to fix M±. The masses of the remaining three 
neutralinos receive one-loop quantum corrections. In this scheme, these counterterms are [28] 

1 / 5m-+ 5m-+ 

bM i = TT2 2 ( M 2 m l+ - »detX)^±- + {M 2 m\ + - fxdetX)- ' 

M 2 — fJ, \ x i rn~+ x 2 



m xt 



L w 

1 / 5m-+ 5m~+ 

S V = ~2 TUt. (» m l+ ~ M 2 detX)^±- + (jj,ml+ - M 2 detX)^±- 

H 2 - Mi \ Xi m-+ x 2 m-+ 

r2 



- M^(/i + M 2 s 2 p) 6 -^f - M 2 M^s 2p c 2 ^) , (3) 

M w H / 

*Mi = ■T=2-(<$m x o - iV 2 VM 2 + 2iV 3i iV4^/i 

-2^1^3*^13 - 2N 2 iN 3i 5Y 2 3 ~ 2N li N 4i 5Y 14 - 2N 2i N M 5Y 24 ) , (4) 

with detX = M 2 /i — M^s 2 /s [28] . <5m^o is the counterterm of zth neutralino defined entirely 
from its self-energy, see[25]. 50 represents the shift on the parameter O that generates the 
counterterm for that quantity. 

Looking at these equations some remarks can be made. First, in the special configuration 
M 2 ~ ±/i an apparent singularity might arise. Ref. [10] pinpointed this configuration which can 
induce a large ^-scheme dependence in the counterterms SMi t2 and 5fJL. Such mixed scenario is 
not covered here. 

Second, the choice of m^o as an input parameter is appropriate only if the lightest neutralino 
is mostly bino (|iV"ix| ~ 1) or if the bino like neutralino is not too heavy compared to other 
neutralinos. Indeed we can see that if Nu ~ the counterterm 5M\ is subject to large un- 
certainty and may introduce large finite correction, this is related to the fact that M\ is badly 
reconstructed. To avoid such uncertainty we only choose i as the most bino like, in other words 
in Eq.gl |iVH| =Max(\N lj \),j = 1, ..4. 

v) Finally diagonal field renormalisation is fixed by demanding that the residue at the pole of 
the propagator of all physical particles to be unity, and the non-diagonal part by demanding 
no-mixing between the different physical particles when on shell. This is implemented in all the 
sectors. In our case apart from the muon, this step is important for the Xi- We insist that Nij 
is used both at the tree-level and one-loop level. Nonetheless to define the physical state we do 
introduce the shift for the neutralinos [28] through wave function renormalisation 



X°^X° + ~£ ( 6Zi i PL + 5Z h P R) X°j • (5) 
3 

vi) Dimensional reduction is used as implemented in the FFL bundle at one-loop through the 
equivalent constrained dimensional renormalisation [33] . 

2.2.3 Infrared divergences 

For the processes X1X1 ~~ * ; we can decompose the one- loop amplitudes in a virtual part 

JWSL and a counter-term contribution Mct- The sum of these two amplitudes must be 
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ultraviolet finite and gauge independent. Due to the virtual exchange of the massless photon, 
this sum can contain infrared divergencies. This is cured by adding a small mass to the photon 
and/or gluon, A 7 and X g . This mass regulator should exactly cancel against the one present 
in the final state radiation of a photon. The QED contribution is therefore split into two 
parts : a soft one where the photon energy E~ is integrated to less than some small cut-off 
k c and a hard part with E 7 > k c . The former requires a photon mass regulator. Finally the 
sum Mii 00 p + Mct + M s ° it {E 1 < k c ) + M^ a g d (E^^ g > k c ) should be ultraviolet finite, gauge 
invariant, not depend on the mass regulator and on the cut k c . 

2.2.4 Checking the result 

i) For each process and set of parameters, we first check the ultraviolet finiteness of the results. 
This test applies to the whole set of virtual one-loop diagrams. The ultraviolet finiteness test is 
performed by varying the ultraviolet parameter Cjjy = 1/e, e is the usual regulator in dimen- 
sional reduction. We vary C\jy by seven orders of magnitude with no change in the result. We 
content ourselves with double precision. 

ii) The test on the infrared finiteness is performed by including both the loop and the soft 
bremsstrahlung contributions and checking that there is no dependence on the fictitious photon 
mass A 7 . 

iii) Gauge parameter independence of the results is essential. It is performed through a set of 
the eight gauge fixing parameters based on the implementation of a non-linear gauge [27] . 

2.3 Effective couplings for neutralino interactions vs Full calculation 
2.3.1 Contributions at full one-loop 




(a) (b) (c) 



Figure 2: Different types of corrections appearing at one-loop for the process XiXi ~~ ^ ^ + 1 1 - ( a ) 
is the correction to vertex f or the t-channel fx exchange, (b) is the full correction to the 

XiXiZ vertex for the s-channel Z exchange and (c) is an example of a box loop. 

The full set of one-loop contributions to the process XiXi ~~ * consist of two-point 

functions (self-energies and transitions such as Xi ~^ X%)i vertex three-point functions as in 
Figs[2^a,b) and box diagrams. The vertex corrections include also the counterterms, the latter 
as explained previously involve two-point functions. To these, one should also add the QED 
final state radiation. 
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2.3.2 Effective couplings of the neutralino at one-loop 

Among this full set of corrections one can construct a finite subset that is not specific to the 
muon. This subset will be involved in all processes involving neutralinos. For example, the ver- 
tex correction to XiXi^ is obviously independent of the muon being in the final state, a similar 
statement can be said for XiXi^/H/A. Also, all occurrences of the wave function renormalisation 
of the neutralino (including transitions between neutralinos) and the Z are process independent. 
The same can be said also of the counterterms to the gauge couplings and the vacuum expec- 
tation values or in other words v,tp. On the other hand the wave function renormalisation of 
the muon is specific to the muon final state. The boxes the four-point one particle irreducible 
(1-PI) functions, as well as the QED correction are also specific to the process. The construct 
of the universal correction for the effective coupling Xiff from Xi^f 1 is different from that of 
XiXi^i since in the latter all three particles can be considered as universal. For example the 
full correction to the vertex Xi^fi shown in Fig. EJa), consists of a 1-PI 3-point function vertex 
correction (triangle) which is muon specific and that does not need to be calculated to build up 
the effective coupling. It also contains wave function renormalisation of the neutralinos as well 
as counterterms for the gauge couplings and for other universal quantities such as tg which must 
be combined to arrive at the universal correction for the Xiff vertex. The aim of the paper 
is therefore to extract these process independent contributions and define effective vertices for 
the LSP interactions. This is akin to the effective coupling of the Z to fermions where universal 
corrections are defined. Describing the bulk of the radiative corrections in terms of effective 
couplings has been quite successful to describe for example the observables at the Z peak. Al- 
though not describing most perfectly the effect of the full corrections for all observables (for 
example Zbb receives an important triangle contribution due to the large top Yukawa coupling) 
one must admit that the approach has done quite a good job. Most of the effective corrections 
were universal, described in terms of a small set of two-point functions of the gauge bosons. 
The other benefit was that such approximations were sensitive to non decoupling effects that 
probe higher scales (top mass and Higgs mass). The set of two-point functions, and for XiXi^ 
three-point functions, should of course lead to a finite and gauge invariant quantity. Loops 
involving gauge bosons have always been problematic (even in the case of the Zff) in such an 
approach since it is difficult to extract a gauge independent value. For the couplings of the neu- 
tralinos as would be needed for approximating their annihilation cross section independently of 
the final state, one would therefore expect that apart from the rescaling of the gauge couplings 
which can be considered as an overall constant, the mixing effect between the different neutrali- 
nos should be affected. One can in fact re-organise a few of the two point functions (that can 
be written also as counterterms) to define an effective coupling for the neutralino. One should 
of course also correct in this manner the coupling. Let us stress again that in this first 

investigation we will primarily take into account the effects of fermions and sfermions in the uni- 
versal loops. For the XiXi^ effective we also attempt to include the virtual contribution of the 
gauge bosons especially that for the higgsino-like the coupling to W and Z are not suppressed. 

2.3.3 The effective Xiff 

To find the process independent corrections to this coupling, we recall that in the basis (B°, W°, H®, H®) 
before mixing and for both fi,R the couplings for the two chiral Lorentz structures writes as 

1 / , v 3 \ 1 i i v 3 9m u gm d , , , g g , 

T2 { 9 Y f , gr f , nj , y 2J ) = -j= (g Y, , gr f , — , _) ^( g>g ,—,—), (6) 

Yf,r£ are the isospin and SU(2) charges of the corresponding fermion/sfermions. The two 
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higgsinos couple differently to the up and down fermions with a coupling that is proportional 
to the Yukawa coupling. Though this is not universal we can still isolate a universal part 
where there is no reference to the final fermion/sfermion. This is what is meant by the last 
expression in Eq. [6] where the explicit mass of the corresponding fermion masses has been 
dropped. The variations/counterterms on these parameters have to be implemented before 
turning to the physical basis. The latter as explained in the previous paragraph is achieved 
through the diagonalising matrix N (Eq. [2|) as in tree-level supplemented by wave function 
renormalisation which involve both diagonal and non diagonal transitions of the neutralino, see 
Eq. [5j In the case of effective coupling of neutralinos, this is achieved by defining an effective 
mixing matrix such that N — > N + AN X ? * in all couplings of the neutralino. The AN X ? * write 
as 



9 1 j 

9 1 j 

where j runs from 1 to 4 and for the LSP , i = 1. 

All the counterterms above are calculated from self-energy two-point functions and are fully 
defined in [271 ES]- Sg/g = Se/e — Ssw/sw^g' /g' = Se/e — Scw/cw- Ssp/sp = ciStp/tp. Eq.0 
agrees with what was suggested in [33]. Let us stress again that in these self-energies no gauge 
bosons and therefore no neutralinos and charginos are taken into account but just sfermions 
and fermions, otherwise this would not be finite. For a bino-like, self-energies containing gauge 
and Higgs bosons (with their supersymmetric conterparts) are not expected to contribute much. 
This is not necessarily the case for winos and higgsinos. 



2.3.4 The effective jfijfiZ 

Since all particles making this vertex can now be considered as being process independent 
(as far as neutralino annihilations are concerned), all counterterms including wave function 
renormalisation of both the Z and Xi must be considered. The price to pay now is that the 
genuine triangle vertex corrections XiXi^ must also be included. It is only the sum of the 
vertex and the self-energies that renders a finite result. When correcting this vertex one must also 
correct the Z(i + fi~ vertex keeping within the spirit of calculating the universal corrections. This 
can be implemented solely through self-energy corrections (excluding the muon self-energies) and 
there is no need to calculate here the genuine vertex corrections. An exception would be the 
production of the b and to some extent the top where genuine vertex corrections are important. 
Talking of heavy flavours, when computing the correction to the XiXi^ with the Z off shell 
with an invariant mass Q 2 , one should also include the X1X1G vertex, where G° is the neutral 
Goldstone boson. In our case we restrict ourselves to almost massless fermions. The case of 
the top and bottom final states will be addressed elsewhere together with the potential relevant 
contribution of the Higgses in the s-channel. 

Since one is including the genuine 1-PI vertex correction, it is important to inquire whether 
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this correction generates a new Lorentz structure beyond the one found at tree-level. The 
contribution to the tree-level Lorentz structure is finite after adding the self-energies and the 
vertex. Any new Lorentz structure will on the other hand be finite on its own. General arguments 
based on the Major ana nature of the neutralinos backed by our numerical studies show that no 
new Lorentz structure is generated for neutralinos. First of all, at tree-level one has only one 
structure 

£xWz = ^(Ni 3 N 13 -N m N u )x iW X°iZ», gz = ^—. (8) 

The overall strength is a consequence of the fact that the coupling emerges solely from the hig- 
gsino with a gauge coupling. Indeed in the (B°, W°, H±, H®) basis the coupling is oc gz(0, 0, 1, — 1). 
Only the Lorentz structure 7^75 survives as a consequence of the Majorana nature. With pi,P2 
denoting the incoming momenta of the two Xii at one- loop a contribution (p^ — p%) does not 
survive symmetrisation, whereas +p%) will not contribute for massless muons. We calculate 
this correction for a Z with an invariant mass Q 2 , in the application this Q 2 will be set to the 
invariant mass of the muon pair. This vertex contribution is denoted Ag^,_ n „(Q 2 ) = Aq^,- n „. 
The contribution of the coupling counterterms defining gz and the Z wave function renormal- 
isation define the universal correction to the Z coupling strength g z = gzO- + &9z), with 
Agz/dz = figz/dz + SZzz/2. 8Zzz is the wave function renormalisation of the Z. We of 
course have to add the wave function renormalisation of the Xi like what was done for the Xiff 
vertex. We improve on this implementation by taking into account the fact that the Z is off- 
shell and therefore the wave function renormalisation through dZzz = N-' ZZ (M Z ) i s onr y part 
of the correction that would emerge from the correction to the complete Z propagator in the 
s-channel contribution with invariant mass Q 2 . Note that here there is no need for including a 
Zj transition since photons do not couple to neutralinos. Collecting all these contributions, the 
effective vertex is obtained by making 

gz -> g% 0z and N iX -> N a + AN§*° Z with 

9xhz = gz(i + Agz(Q 2 ) + Ag* Mz (Q 2 )); (9) 
AN mz = lj2N kj SZ ki , (i,j,k) = l..A. (10) 



k 



Explicitly 



1 f-n' m\ o s ^n 7Z (0)\ , 1/ c 2 w \ fU Z z(M 2 ) TlwwiM 2 ,] 



*9Z = ^ n' (0)-2^^p + 



2 V 77 c w M 2 J 2 V s 2 w ) V Mf Mfr 

1 ( U ZZ (Q 2 ) - n Z z(M 2 ) 

2 V Q 2 - M 2 



fill 



At the same time for the fermion with charge qj we correct the Zff ex <?z(75 + (1 — ^\Qf\ s w))liJL 
by effectively making gz — > gz(l + Agz) with Agz defined in Eq. [11] and s 2 ^ to 

As 2 = fjV ( n Z z(M 2 ) _ U WW (M^) \ + c w U yZ (k 2 ) (12) 



,W s 2 w \ Mf ~ M 2 W J^s w k 2 

By default we include only the fermions and sfermions in the virtual corrections described by 
Eqs. [Tlll2i For the XxXx^ one ex P ec t the contribution of the gauge bosons and the neutrali- 
nos/charginos to be non negligible especially for the higgsino case. In fact, including such 
contributions still gives an ultraviolet finite result for g_^_ „ in Eq. [11] which is a non trivial 
result. Moreover this contribution is gauge parameter independent in the class of (linear) and 
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non- linear gauge fixing conditions [29 . To weigh up the gauge/gaugino/higgsino contribution 
we will therefore also compare with this generalised effective g e ^- 0z including all virtual parti- 
cles. Observe that in Eq. [11] we have the contribution TLyz(0) which vanishes for fermions and 
sfermions but which is essential for the contribution of the virtual W. In any case including 
gauge bosons in the renormalisation of electromagnetic coupling requires the inclusion of the 
n 7 z(0) in Eq. [TT]for gauge invariance to be maintained [29]. We stress that we will present the 
effect of the generalised effective coupling g e ^ _ „ as an indication of the gauge boson contribution 
while keeping in mind that this result may lead to unitarity violation. Indeed through cutting 
rules, the W loop can be seen as made up of the scattering W + W~ — > Z — > [i> + fJ>~ that needs a 
compensation from the cut in the box shown in Fig. [2jc). For the effective Z[i + fi~ coupling we 
only include the fermion/sfermion contribution in Eqs. [1111121 adding the gauge bosons would 
require part of the 1-PI triangle contribution to Z — > n + /i~ . 



3 Analysis 

Since we will be studying different compositions of the neutralinos we will take different values 
for the set Mi,M%, fj,. On the other hand the default parameters in the Higgs sector are 

M A0 = ITeV tp = 4. (13) 

The sfermion sector is specified by a rather heavy spectrum (in particular within the limits 
set by the LHC for squarks[lj). All sleptons left and right of all generations have a common 
mass which we take to be different from the common mass in the squark sector. All tri-linear 
parameters Af (including those for stops and sbottom) are set to 0. The default values for the 
sfermion masses are 

M Ir = M Jl = 500GeV, 
Mul = M~ r = Mq l = 800GeV, 
A f = 0. (14) 

By default we will focus on relatively light neutralinos (around 100 GeV) scattering with a 
relative velocity v = 0.2. 

To analyse consistently the efficiency of effective corrections we will refer to the following 
quantities : 



A e ff = • (15) 

Here a e g is the cross section calculated with the effective couplings that include, by default, 
universal process independent particles excluding gauge bosons and gauginos/higgsinos. We 
will explicitely specify when including all virtual particles in those corrections, referring to 
it as A^. This correction will be compared to the correction solely due to the running of 
the electromagnetic coupling, see Eq. [TJ To see how well the correction through the effective 
couplings Xiff an d XiXi^ reproduces the full one-loop correction we introduce 

A ^one-loop &c5 
NU = , 

0*0 

A0one-loop 00 
full = • (16) 

00 
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with <7one-ioop the full one-loop cross section, Ajvx/ measures what we will refer to as the non- 
universal corrections although strictly speaking this measures the remainder of all the corrections 
that are not taken into account by the effective vertices approach. Af u u = A e g + Ajvc/ is the 
full one-loop correction. 

3.1 Bino Case 

3.1.1 Effective vs full corrections 

We first take (Mi,M2,/i) = (90,200,-600) GeV which yields a lightest bino-like neutralino 
(the bino composition is 99%) with mass m^o = 91 GeV. At tree-level the cross section for 

relative velocity v = 0.2 is cr ^+ /J - = 6.75 x 10 _3 pb. Note for further reference that this is an 

order of magnitude larger than annihilation into a pair of W's: a^y +w _ = 4.51 x 10~ 4 pb. The 
annihilation proceeds predominantly through the t-channel, binos coupling to Z are very much 
reduced. This leads to the following set of corrections 

A eff = 17.52%(A a = 14.56%) A NU = 2.06%(A full = 19.58%). (17) 

For our first try the effective universal coupling does remarkably well falling short of only 2% 
correction compared to the full calculation. Note that although the most naive implementation 
through a running of the electromagnetic coupling fares also quite well it is nonetheless 5% off 
the total correction, therefore the effective correction through the effective couplings performs 
better. It must be admitted though that the bulk of the correction is through the running of a. 
To see how general this conclusion is we scanned over the set (Mi, M%, fJ-) while maintaining Xi 
with a 99% bino like component. This is simply obtained by taking M2 = 500, \i = —600 GeV 
and scanning up to M\ = 350 GeV . We also checked how sensitive our conclusion is depending 
on tp by varying tp from 2 to 40. The suspersymmtery breaking sfermion masses were first 
left to their default values. As Fig. [3] shows, our conclusions remain quantitatively unchanged. 
There is no appreciable dependence in tp, we arrive at the same numbers as our default tp 
value. As for the dependence in M\ it is very slight, for M\ ~ 50 GeV there is perfect matching 
with our effective coupling implementation, then as M\ increases to 350 GeV, the non universal 
corrections remain negligible, below 2%. 

The annihilation of neutralinos and hence the relic density is a very good example of the 
non decoupling effects of very heavy sparticles, a remnant of supersymmetry breaking. The 
variation in the fermion/sfermion masses is all contained in the effective couplings that we 
have introduced. Leaving the dependence on the smuon mass at tree-level, and the very small 
(see below) contribution of the smuon to the 1-PI vertex Xil^fi^ the bulk of the smuon mass 
dependence is within the effective coupling. Fig. [3] shows how the correction increases as the 
mass of the squarks increases from 400 GeV to 3 TeV, we take here a common mass for the 
supersymmetry breaking squark masses (both right and left in all three generations). The non 
universal correction of about 2% is insensitive to this change in squark masses whereas both 
A e g and A]\uj show the same logarithm growth that brings a 3% change as the squark mass is 
varied in the range 400 GeV to 3 TeV. This result also confirms that genuine vertex corrections 
and box corrections are very small. 

We have also extracted the individual contribution of each species of fermions to the total 
non-decoupling effect of sfermions. To achieve this we numerically extracted the logarithm 
dependence of the non decoupling effect for each species of sfermions. We have parameterised 
the effective correction as 

A? = ajlnnij/Q — aflnrrif/Q + bf with Q = 2m^o / = d,R + ur-\- Ql (18) 
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Figure 3: Corrections to the tree-level cross-section for the process XiXi ~ > * n the bino 

case as a function of M\ (left panel) and tp right panel. We show the full one-loop, the effective 
correction and the difference which we term non effective. M2 = 500, /i = —600 GeV. 
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Figure 4: Corrections to the tree-level cross-section for the process X1X1 ~^ ^ ^~ ^ n the bino 
case (Mi = 90, M2 = 200, = —600 GeV) as function of the common soft supersymmetry 
breaking squark mass. 
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The coefficients of the fit are given in Table [TJ As expected the fit to aj is extremely well 
reproduced by the running of a, i.e, aj = N c qj^. We also find ag = (if = a^,b e = = b T . 
The fit to df is made to validate the fit procedure. The most important observation is that 





Qr, 




d R 


a f 


bf 


e 


0.0010 




0.00231 


0.00310 


0.15% 


CM) 


0.000575 


0.00236 


0.000698 


(0.00413,0.00103) 


0.15% 


(t,b) 


-0.00406 


0.00838 


0.000661 


(0.00413,0.00103) 


0.16% 



Table 1: Coefficients of the ln(mj) (running couplings) \n(mj)(non decoupling effects) in A e g. 
(c,s) give very similar results to (u,d). 

the stops behave differently, this is due to the Yukawa coupling of the top and mixing. If there 
were not a compensation between left and right contribution of the stops (compare to u) the 
contribution of the stops would be even more important and would dominate. Considering the 
different contributions and the scales that enter our calculations it is difficult to attempt at 
giving an analytical result, but leaving the stop aside the different contributions to aj can be 
roughly approximated by y^NcNd/8/cyy, Nd = 2 for doublets and 1 for singlet of SU(2). yf is 
the hypercharge, corresponding to the couplings of the sfermions to the bino component. 

3.1.2 Scheme dependence in the bino case 

We have compared the full correction to an approximate effective implementation and observed 
that the approximation is quite good. However, even the full correction, being computed at 
one-loop, it is potentially dependent on the renormalisation scheme chosen. As discussed earlier 
we analyse the tg scheme dependence and the M\ scheme dependence. For tg we obtain the 
following corrections: 

19.58%(DCPR), 19.79%(0Di2), 19.51%(MF). 

This confirms that the tg scheme dependence is very negligible. For the bino case it is natural 
to reconstruct M\ from the LSP, nonetheless analysing the M\ scheme dependence one chooses 
another neutralino, say x[> which in our example is a wino-like. This introduces more uncertainty 
or error since with this scheme the corrections attain 24.08%, more than 4% compared to the 
usual scheme. 

3.2 Higgsino Case 

3.2.1 Effective versus full corrections 

In the bino case our trial point had a neutralino of mass 91 GeV. We therefore take the point 
(600,500,-100) which gives a LSP with M^o = 95 GeV with a 99% higgsino content. The 
sfermion parameters are the default values. In the higgsino case the cross section is dominated 
by the exchange of the Z in the s-channel, so the bulk of the corrections through the effective 
couplings will be through the effective XiXi%- For further reference note that the tree-level 
cross section for annihilation into muons is cr ^+ /J - = 2.58 x 10~ 3 pb, tiny and totally insignificant 
especially compared to annihilation into W, a^ v+w _ = 18.83 pb. This is an observation we will 
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keep in mind. The one-loop corrections we find for o~^ +fi _ are 
(for (jl = -lOOGeV) A cfr = 13.55(A Q = 14.62%) A NU = -21.09%(A /tlH = -7.54%) (19) 

This result is in a quite striking contrast to the bino case. The effective coupling does not 
reproduce at all the full correction and is off by as much as 21%. It looks like, at least for this 
particular choice of parameters, that going through the trouble of implementing the effective 
XiXi^ was m vain since this correction is, within a per-cent, reproduced by the naive running 
of a. As we will see both these conclusions depend much on the parameters of the higgsino and 
even the squark masses. For example consider /i = —50 GeV, leaving all other parameters the 
same. Of course this is a purely academic exercise, since in this case, the charginos with mass 
m ± = 55 GeV are ruled out by LEP data. Nonetheless, in this case 

(for /x = -50GeV) A cff = 10.7(A Q = 12%) A NU = -6.9%(A /uH = 3.8%). (20) 

Had we included all particles in the effective vertex, we would get a correction A^ = 4.4% 
improving thus the agreement with the one-loop correction for this particular value of /i up to 
0.6%. At the same time a correction in terms of a running of a will be off by more than 8%. 
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li (GeV) 

Figure 5: Corrections to the tree-level cross-section for the process XiXi ~~ ^ * n the hig- 

gsino case as a function of /j,. Shown are the effective vertex correction (Effective, with only 
fermions/ sfermions in the loops), the effective X\X\^ coupling including all particles denoted 
Effective (All), the non QED boxes (Boxes) and the full correction. M2 = 500, /x = —600 GeV. 



These two examples show that one can not, in the higgsino case, draw a general conclusion 
on the efficiency of the effective coupling as what was done in the bino case. Let us therefore 
look at how the corrections change with /i, and therefore with the mass of the LSP, while 
maintaining its higgsino nature. We have varied \x from —200 GeV to —40 GeV . Fig. [5] shows 
that the full correction is extremely sensitive to the value of \x. For \x = —200 GeV the full one- 
loop correction is much as —42%, casting doubt on the loop expansion. The effective coupling 
corrections with only fermions/sfermions on the other hand is much smoother and positive 
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bringing about 10% correction. Including all particles in the effective XiXi^ vertex brings in 
an almost constant reduction of about 6%. Therefore as the value of increases the effective 
one-loop corrections in the case of the higgsino case can not be trusted. The same figure shows 
that the behaviour and the increase in the corrections is due essentially to the contribution of 
the boxes. Here the boxes mean the non QED box (involving an exchange of a photon which 
are infrared divergent before including the real photon emissiorJl) . The large contribution of the 
boxes can be understood by looking at the box in Fig. [2(c). Indeed, as argued previously, cutting 
through the box reveals that it represents x?X? ~~ > W + W~ production that rescatter into y + y~ ■ 
Both these process have a very large cross sections compared to the tree-level X1X1 ~^ / U+ M - 
Our conclusion is therefore that the effective vertex approximation is inadequate as soon as 
the channel XiXi ~~ * W + W~ opens up. When this occurs, in practical calculations of the relic 
density, the channel x?X? ~~ * is irrelevant and must rather analyse the loop corrections to 

XlXi ~^ W + W~. This process was studied in [131 EH] and will be investigated further through 
an effective approximation in a forthcoming study. 

On the other hand, the dependence of the relative correction on tp is quite modest even though 
there is certainly more dependence than in the bino case, especially at lower values of tan/3. 
This is shown in Fig. [6[ 



Full One-Loop 
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Non-Effective 
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Figure 6: Corrections to the tree-level cross-section for the process XiXi ~~ ^ * n ^ e higgsino 

case as a function of tp. We show the full one-loop, the effective correction and the remainder 
(Non-effective), y. = -100, M 2 = 500,^ = -600 GeV . 



We now investigate the non-decoupling of very heavy squarks (and heavy sfermions in gen- 
eral). Since we are in a Higgsino scenario we expect the Yukawa of the fermions to play a more 
prominent role than what was observed in the bino case. This is well supported by our study. 
Fig. [7] shows how the effective (with only fermions and sfermions) and the full correction gets 
modified when the common mass of all squarks (all generations, left and right) increases from 
400 GeV to 3 TeV. To better illustrate the important effect of the Yukawa of the top/stop sector 
we plot the corrections also for m± = 0.1 GeV. For mt = 170.9 GeV, the correction drops by 
about 13% when the mass of the squarks increase from 400 GeV to 3 TeV. This is much more 

2 The contribution of the QED box + real photon emission is only 0.1% 
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dramatic than in the bino case where we observed a 3% increase in the same range. Observe that 
for our default squark mass of 800 GeV, the effective correction including sfermions/fermions is 
such that it almost accidently coincides with the running of a. If one switches off the top quark 
mass, instead of a 13% decrease we observe an 8% increase for m t = 0.1 GeVl Observe that 
the difference one sees for niQ = 400 GeV between nit = 170.9 GeV and mt = 0.1 GeV is due 
essentially to the running of a with very light top that accounts for 3%. 




400 800 2000 3000 400 800 2000 3000 

Mq (GeV) Mq (GeV) 



Figure 7: Corrections to the tree-level cross-section for the process XiXi ~^ m the higgsino 

case (Mi = 600, M2 = 500, /i = —100) as function of the common squark mass. The right panel 
illustrates the case mt = 0.1 GeV . 

The special role played by the top can be seen even more clearly from each individual 
contribution of the fermion/sfermions and the fit of the contribution according to Eq. [18] as was 
done for the bino case. The contribution of the stop is clearly (especially through Q£) an order 







a>u R 


dft 


a f 


*/ 


e 


0.00304 




0.000366 


0.00309 


-0.12% 


(u,d) 


0.00861 


0.000489 


0.000122 


(0.00414,0.00101) 


-0.13% 


(t,b) 


-0.0701 


0.000826 


0.000108 


(0.00414,0.00101) 


0.13% 



Table 2: Coefficients of the ln(m/) (running couplings) \n(mj)(non decoupling effects) in A e g. 
(c, s) give very similar results to (u,d). Higgsino case. 

of magnitude larger than for all other sfermions, see Table El It is the only one that brings 
a negative contribution. Since this effect is in the universal XiXi^ h wu l show up in many 
processes where the higgsino contributes. 

3.2.2 Scheme dependence in the higgsino case 

We analyse here the tp scheme dependence and the M\ scheme dependence. For tp we obtain 
the following corrections: 

-7.5%(DCPR), -12A%(0DR), -4.76%(MH). 

As expected and in line with the behaviour of the corrections with respect to tan /3, Fig. [6l we see 
that the corrections though larger than in the bino case are nonetheless within 5%. On the other 
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hand, expectedly the choice of Mi has less impact than in the bino case where the reconstruction 
of Mi is essential to define the LSP. In the case of the higgsino, changing the Mi scheme turns 
the full correction from -7.5% (in DCPR scheme for tp) to -10.7%, a 3% uncertainty. 

4 Conclusions 

Very few analyses have been done taking into account the full one-loop corrections to the an- 
nihilation cross sections entering the computation of the relic density despite the fact that this 
observable is now measured within 3% precision. In supersymmetry radiative corrections have 
been known to be important, yet practically all analyses that constrain the parameter space of 
supersymmetry are performed with tree-level annihilation cross sections. Taking into account 
the full one-loop corrections to a plethora of processes is most probably unrealistic. On the 
other hand one must incorporate, if possible simply and quickly, a parameterisation of the the- 
ory error or implement the corrections through effective couplings of the neutralino, in the case 
of supersymmetry. This is what we have attempted in this study for two of the most important 
couplings of the neutralinos Xiff an d XiX?^- I n order to look more precisely at the impact of 
each of these effective couplings we take as a testing ground a most simple process, XiXi ~~ ^ 
and select a neutralino that is either almost pure bino or pure higgsino. We do not strive at 
finding a scenario with the correct relic density since our primary task is to study this vertices 
and the approximations in detail. In this exploratory study taking a final state involving gauge 
bosons would only confuse the issues. Nonetheless, the impact of the gauge bosons is studied. 
Indeed, we have shown how the construction of the effective XiXi^ is quite different from that 
of the Xiff- For the latter the effective coupling involves self-energy corrections, whereas for 
the former the one-particle irreducible vertex correction must be added. These examples and 
the construction of the effective coupling already pave the way to a generalisation to the effec- 
tive couplings XiXih, H, A and XiX + W which we will address in forthcoming publications with 
applications to different process, including gauge boson final states. Even with the effective 
couplings we have derived, we could generalise the study of XiXi ~~ ^ to cover not only 

pure winos, but also mixed scenarios and also heavy fermions. 

Our preliminary study on the simple process XiXi ~* ^ >~ is already very instructive. To 
summarise the bino case, we can state that the effective couplings approach is a very good 
approximation that embodies extremely well the non decoupling effects from heavy sfermions, 
irrespective of many of the parameters that are involved in the calculation, as long as one is 
in an almost pure bino case. The effective coupling implementation is within 2% of the full 
one-loop calculation. Here, this reflects essentially the correction to the Xiff coupling. The 
scheme dependence from tp is very small, this result stands for large Mi masses as long as the 
neutralino is more than 90% bino like. In particular for higgsino-like LSP in excess of 90 GeV as 
imposed by present limits on the chargino, the effective coupling implementation in the annihi- 
lation XiXi ~~ * /• l+ A t- fails. It worsens as the mass increases due to the importance of a large box 
contribution corresponding to the opening up of XiXi ~~ ^ W + W~ which would in any case be 
the dominant process to take into account when calculating the relic density. The large Yukawa 
of the top has a big impact on the radiative corrections and in particular on the non-decoupling 
contribution of a very heavy stop. Although this is an example which shows, in principle, the 
failure of the effective approach apart from correctly reproducing the non-decoupling effect of 
very heavy squarks, we need further investigation on the dominant processes, in this case anni- 
hilations into W, Z, to see if these dominant processes could on the other hand be reproduced 
by an effective coupling approach. If the effective approach turns out to be efficient for the 
dominant processes, where and if the box corrections are tamed, the effective coupling could 
still be a good alternative for the calculation of the relic density with high precision. We leave 
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many of these interesting issues to further analyses. 
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